home *** CD-ROM | disk | FTP | other *** search
/ IRIX 6.2 Development Libraries / SGI IRIX 6.2 Development Libraries.iso / dist / complib.idb / usr / share / catman / p_man / cat3 / complib / cunghr.z / cunghr
Text File  |  1996-03-14  |  2KB  |  67 lines
  1.  
  2.  
  3.  
  4. CCCCUUUUNNNNGGGGHHHHRRRR((((3333FFFF))))                                                          CCCCUUUUNNNNGGGGHHHHRRRR((((3333FFFF))))
  5.  
  6.  
  7.  
  8. NNNNAAAAMMMMEEEE
  9.      CUNGHR - generate a complex unitary matrix Q which is defined as the
  10.      product of IHI-ILO elementary reflectors of order N, as returned by
  11.      CGEHRD
  12.  
  13. SSSSYYYYNNNNOOOOPPPPSSSSIIIISSSS
  14.      SUBROUTINE CUNGHR( N, ILO, IHI, A, LDA, TAU, WORK, LWORK, INFO )
  15.  
  16.          INTEGER        IHI, ILO, INFO, LDA, LWORK, N
  17.  
  18.          COMPLEX        A( LDA, * ), TAU( * ), WORK( LWORK )
  19.  
  20. PPPPUUUURRRRPPPPOOOOSSSSEEEE
  21.      CUNGHR generates a complex unitary matrix Q which is defined as the
  22.      product of IHI-ILO elementary reflectors of order N, as returned by
  23.      CGEHRD:
  24.  
  25.      Q = H(ilo) H(ilo+1) . . . H(ihi-1).
  26.  
  27.  
  28. AAAARRRRGGGGUUUUMMMMEEEENNNNTTTTSSSS
  29.      N       (input) INTEGER
  30.              The order of the matrix Q. N >= 0.
  31.  
  32.      ILO     (input) INTEGER
  33.              IHI     (input) INTEGER ILO and IHI must have the same values as
  34.              in the previous call of CGEHRD. Q is equal to the unit matrix
  35.              except in the submatrix Q(ilo+1:ihi,ilo+1:ihi).  1 <= ILO <= IHI
  36.              <= N, if N > 0; ILO=1 and IHI=0, if N=0.
  37.  
  38.      A       (input/output) COMPLEX array, dimension (LDA,N)
  39.              On entry, the vectors which define the elementary reflectors, as
  40.              returned by CGEHRD.  On exit, the N-by-N unitary matrix Q.
  41.  
  42.      LDA     (input) INTEGER
  43.              The leading dimension of the array A. LDA >= max(1,N).
  44.  
  45.      TAU     (input) COMPLEX array, dimension (N-1)
  46.              TAU(i) must contain the scalar factor of the elementary reflector
  47.              H(i), as returned by CGEHRD.
  48.  
  49.      WORK    (workspace/output) COMPLEX array, dimension (LWORK)
  50.              On exit, if INFO = 0, WORK(1) returns the optimal LWORK.
  51.  
  52.      LWORK   (input) INTEGER
  53.              The dimension of the array WORK. LWORK >= IHI-ILO.  For optimum
  54.              performance LWORK >= (IHI-ILO)*NB, where NB is the optimal
  55.              blocksize.
  56.  
  57.      INFO    (output) INTEGER
  58.              = 0:  successful exit
  59.              < 0:  if INFO = -i, the i-th argument had an illegal value
  60.  
  61.  
  62.  
  63.                                                                         PPPPaaaaggggeeee 1111
  64.  
  65.  
  66.  
  67.